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ABSTRACT 



We present Herschel/SPIRE images at 250, 350, and 500 yum of NGC 4244, a typical low-mass, disk-only and edge-on spiral galaxy. 
The dust disk is clumpy and shows signs of truncation at the break radius of the stellar disk. This disk coincides with the densest 
part of the Hi disk. We compare the Spectral Energy Distribution, including the new SPIRE fluxes, to 3D radiative transfer models; 
a smooth model disk and a clumpy model with embedded heating. Each model requires a very high value for the dust scale-length 
{hj = 2 - 5 h,), higher dust masses than previous models of NGC 4244 (Mj = 0.47 - 1.39 x 10^ Mq) and a face-on optical depth 
of T^y"' = 0.4 - 1.12, in agreement with previous disk opacity studies. The vertical scales of stars and dust are similar The clumpy 
model much better mimics the general morphology in the sub-mm images and the general SED. The inferred gas-to-dust mass ratio is 
compatible with those of similar low-mass disks. The relatively large radial scale-length of the dust disk points to radial mixing of the 
dusty ISM within the stellar disk. The large vertical dust scale and the clumpy dust distribution of our SED model are both consistent 
with a scenario in which the vertical structure of the ISM is dictated by the balance of turbulence and self-gravity. 

Key words. Radiative transfer, (ISM:) dust, extinction, ISM: structure. Galaxies: ISM, Galaxies: spiral. Galaxies: structure 



1. Introduction 

The edge-on perspective of a spiral galaxy's disk reveals both 
the vertical structure of the disk and its faint outskirts, thanks 
to line-of-sight integration. One typical feature is the mid-plane 
dust absorption lane. Dust is linked to the cold molecular ISM, 
mechanically (Allen et al., 1986; Weingartner & Draine, 2001), 
through mutual s hielding, and a s a catalyst for molecular gas 
formation. .Dalcanton et al.l (12004) used the appearance of dust 
lanes as a probe of vertical stability of spiral disks. They found 
that in massive spirals, the ISM collapses into a thin dust-lane, 
while in low-mass disks (v,ot < 120 km/s) the dust morphology 
is flocculant. 

One comprehensive approach to investigate this phe- 
nomenon is to model the multi-wavelength information on a 
range of spiral disks with a Spectral Energy Distribution (SED) 
model. Dust absorption, especially within the thin dark dust 
lane, is strongest at optical and ultra-violet wavelengths and 
the challenge for the models is to balance the observed ab- 
sorption with the dust emis sion at far-infrared (FIR) and sub 



millimet er waveleng t hs (e.g ., Popescu et al., 2000; Misselt et al 



engi , , 

2001; .Gordon et aU 120011: Alton et aL .2004.; Bianchj, .2008 



Baes et all 120101; iBianchi & Xilourir i201 ih . Until recently. 



these models were degenerate in the vertical distributions of stel- 
lar light and ISM due to a lack of resolution and wavelength 
coverage in the FIR and sub-millimeter, but with the advent of 
the Herschel Space Observatory, the vertical structure of nearby 
edge-on disks can now be fully resolved. 

Massive edge-on spiral ga laxies have been obs e rved by sev 
eral H e rschel survey s (see iBoselli et al.l 120101: lEales et al " 
l20T0bl iDavies et all I20T0I) and a dedicated prograi 
(Verstappen et al. in prep.) which targets seven massive 
spirals. Our NHEMESES0 program aims to observe a comple- 
mentary sample of low-mass nearby edge-ons. The combined 
observations will serve as a test for the different radiative 
transf e r models in the literature; those fro m Baes et al. (I2003L 
20TT), Bianchi' (^2007, '200?), Gordon et al.' ('2001');'Misselt et alj 
(2001) and Pierini et al. (2004), and Popescu et al. (2000, 20 ij^ 
These models will further test the suggested dichotomy in the 
vertical dust stiTicture. 

Here, we explore the first NHEMESES results on NGC 
4244, the prototypical low-mass (Vrot - 95 km/s), late-type 
(Sc), edge-on spiral at a distanc^l of 4.4 Mpc to compare to the 
first Herschel result s on th e massive ed ge-on spiral NGC 891 
(iBianchi & Xilourisl 1201 lb . Oiiginally, Ivan der Kruit & Searld 
(|1981|) found a single edge-on exponential disk truncated at 9 
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' HERschel Observations of Edge-on Spirals, HEROES 

- New HErs chel Multi-wav elength Extragalactic Survey of Edge-on 

Spirals (Holwerda et al.L 1201 ih 

^ The mean of the dista nces in NED from iHeald et^ i201 Ibl) . 

iRadbum-Smith et al derive a distance of 4.7 Mpc. 
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kpc, later confirmed by iFrv et al.l (Il999l) . The t runcation was 
also observed with the Hubble Space Telescope (Ide Jong et al.L 
[2007). Yoachim & Dalcanton (2006) n ote a s econd, thick stel- 
lar disk as doubtful, but lComeron et akl (1201 ih report both a thin 
and thick di sk of equa l mass . Its center hosts a rotating nuclear 
star-cluster ( Seth et al. , l2008b . and the disk is classified as "cor- 



contained within the D25 radius, and the profile shows clumps 
of emission at both ends of the disk. These clumps of sub- 
mm emission are the last reliable det ection of the disk be- 
fore the truncation radius (420", 9 kpc Ivan der Kruit & SearleL 
1981). The vertical sub-mm profiles are similar in width to 
those in Spitz er 3.6 ixm emissio n, a good tracer of the stel- 
lar mass (e.g. iMeidt et all l2012h . Figure |2] suggest that dust 
is confined to the thin disk, in contrast to NGC 891, where 
a second vertical component is observed (Kamphui s et all 
12007'; 'Bianchi & Xilouris', '201 1' ), as well as a thin dust lane 
(Xilouiis et al., 1999; Schechtma ri-Rook et al.Ll2012allbl) . 

. , Fig. |3]shows the optical and infrared SED of the galaxy. The 

IZschaechner et al.L l20lU iHeald et all 12011 ah. The Hj mass , SPIRE flux densities show a spectral index Fy oc v' "', flatter than 

the more massive NGC 891 {Fy oc v^^), and thus indicative of 



rugated' '(Florido et allll991h 



NGC4244's ISM is m ostly in a thin disk, both the ionized 
gas dHoopes et all 1 19991) and warm dust, implying a low spe- 
cific star-formation rate (iKodaira & Yamashitai Il996l) . The Hi 
disk shows a clear warp, kinematic evidence for a compact 
dark matter halo, and la g ging extra-p l anar ga s (lOlUngi 119961: 



based on CO observations, is 1.4 x 10^ Mq dMatthews & Wood 

l200l . 

A recent SED model bv lMacLachlan et al.l (I2OIII) finds that 
most (80%) of the FIR emission could arise from an optically 
thi n disk with a sim ilar scale-height for both dust and stars (as 
did lSeth etalll2005L based on stellar populations), a dust scale- 
length 1.8 times the stellar one and a dust mass of 2.38 x IO^'Mq. 
They note their model diverges from the observed disk in the 
FIR and that the dust model is still uncertain without sub-mm 
observations. 

In this Letter, we present new SPIRE sub-mm observations 
of NGC 4244, adding spatial information on the cold dust, and 
present two TRADING models to illustrate the role of ISM ge- 
ometry on the observed overall SED and the sub-mm images. 

2. Observations 

SPIRE (Spectral and Photometric Imaging REceiver, 
iGriffin et all 1201 Ol) instrument at 250, 350 and 500 u.m 
onboard the ESA Herschel Space Observatory (Pilbra tt et all 
[20 10) observed NGC 4244 at 250, 350 and 500 pm in large 
map mode, over 25'x25' centered on the object with two 
cross-scans, using a 30" s ~' scan rate. Data-reduction is as in 
iBianchi & XilourisI (1201 ll) : we generated maps with the naive 
map-making procedure within HIPE (Ott, 2010), with pixel 
sizes of 6", 8", 12" at 250, 350 and 500 pm, respectively (about 
1/3 of the measured FWHM: 18", 25" and 36"), resulting in 
uniform backgrounds and sky rms noise of 1, 0.6 and 0.3 MJy 
sr~' respectively. Fig. [U. A color correction half-way between 
that for a point- and extended-source was implemented with 
a color-correction in between the point- and extended-source 
(marginally resolved vertically), we derive integrated fluxes 
of 15.9, 9.8, and 5.1 Jy (250, 350 and 500 pm), and adopt a 
conservative 15% calibration error 



3. Results & Analysis 

Figure [1] shows a SPIRE co lor-compo site, together with the Hi 
contours from Zscha echne r et al.l (1201 1) . The dust emission falls 
entirely within the second highest Hi contour, concentrated in a 
flat, single disk with no emission in the extended warped Hi disk 
and envelope, in contra st to the more massiv e NGC 891, where 
iPopescu & Tuffs! (I2003h and lBianchi & XilourisI ( 120111) find evi- 
dence for dust througho ut the Hi disk and the extended envelope 
dOosterloo et alll2007l) . Dust is concentrated in the thin Hi disk, 
but we note that dust and Hi clumps do not coincide. 

Figure |2] shows the radial and v ertical profiles of the sub- 
mm data, SDSs-R data dAihara et all |201 Ij) . and the Spitzer 3.6 
pm image from a related Spitzer prograrry. The sub-mm flux is 



a colder dust temperature (or of different dust properties). When 
compared to the radiative transfer model of MacLachla n et al.l 
(I2OI ih . the observed fluxes at 350 and 500 pm are higher by a 
factor two (Figurel3l lHolwerda et al.Ll20TTl) . 

We model the surface brightness distribution and SED 
of NGC 4244 using the radiative transfer code TRADING 
^Bianchi, 2008). We initialize the stellar dist ribution w i th the 
scale- length and height (/j,,z*) as derived bv iFrv et al.l (1 19991) 
(scaled to 4.4 Mpc), and adopted an intrinsic SED that, after be- 
ing processed by dust, reproduces the observed stellar fluxes in 
the optical and NIR. We truncate the model, for both stars and 
dust, at 420" (9 kpc) since it is to this radius we have reliable 
fluxes for both. 

First, we try a smooth exponential dust disk with the same 
scale-height {zd) as for stars, and modify the scale-length {hd) 
and the central V-band face-on optical depth though the whole 
disk (r^" ) to approximate both the sub-mm surface brig htness 
profile s (Fig. 13 and the SED (Fig. O. As MacLachla n et al.l 
(I2OI 1). we require a dust scale-length larger than the stellar 
one to match the observed major axis profiles. However, con- 
trary to iMacLachlan et al.l (l201 ll) . we need a higher opacity 
(and hence dust mass) to match the peak of the SEEQ Ours 
and MacLachlan et al. (2011) models are compared in Table [1] 
While the smooth model can broadly match the peak of the SED 
and the profiles, it lacks hot dust to explain emission at IQpm. 
And it under-predicts the 500 pm flux, implying a lower dust 
temperatur^ 

Secondly, we consider an exclusively clumpy disk model; an 
exponential distribution of Giant Molecular Clouds (GM C), sim- 
ilar to the clumpy part of the smooth +clumpy model in Bianchil 
(l2008h with Mgmc - lO ' ^-^Mp, gas-to-dust mass ratio - 100, 
Rgmc - 100 pc. Bianchi (2008) find the SED does not depend 
critically on the GMC size distribution. This model matches the 
sub-mm SED well but fails in the FIR (Fig. [3] dashed line). To 
match the 70 pm emission, we heat 60% of the clouds internally 
with a hot young stai'cluster (T = 35000 K, 1.56 x 1O''L0, Fig. 
|3] solid green line), simil ar to the secon d dense dust disk with 
young stars introduced bv lPopescu et al.l(l2000i) to reconcile FIR 
and sub-mm fluxes. We note that PAH fluxes (e.g., Spitzer 8 pm) 
are not yet reproduced by this model but this may be solved with 
a range of cloud sizes and heating source types in a future iter- 
ation. The resulting heated clumpy model requires a large dust 



P.I. R.S. de Jong, Holwerda et al. in preparation. 



Although the model o f IMacLachlan et al.1 d201 ih includes a distri- 
bution of compact sources, with a dust mass about ten times than that 
of their diffuse disk, its contribution to the peak SED, as estimated from 
their Fig. 12, is smaller than 15%. 

* Our assumption is that the dust grain population of NGC 4244 and 
Milky Way are identical but differences in sub-mm flux and optical ex- 
tinction can be due to different grain distributions as well. 
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Fig. 1. A colour com- 
posite of the SPIRE 
maps of NGC 4244, 
250 (red), 350 (green), 
and 500 fim (blue) with 
the Hi contours from 
IZschaechner et al] ( 120111) 
(contours at 6.4 x 10" 
cm"- increasing by fac- 
tors of 2. The SPIRE flux 
stays within the 8. 1 x 10" 
cm"^ contour (second to 
highest). 
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Fig. 2. Major (left) and minor (right) axis pro- 
files, averaged over a strip of 30" (solid lines), 
in the sdss-r (green) and IRAC 3.6 fim (pink, 
offset by a factor 2 for clarity) and the 250 
(blue) and 500 (red) /im SPIRE data, all at na- 
tive resolution. The major axis profile spans 
£>25 = 1616 (RC3). We also show the smooth 
(dot-dash) and clumpy (dashed) TRADING 
models for comparison (convolved with the 
SPIRE PSF s at 250 and 500 jjm ). We assumed 
PA.=47.3° jde Jong etaUl2007l) . 



Table 1. Results from the SED models. 



Parameter 


MacLachlan 


TRADING 




Unit 




etal. 2011 


Smooth 


Clumped 




h. 


1.9 


2.3 


2.3 


kpc 


Z.,Zd 


0.2 


0.3 


0.3 


kpc 


hd 


3.4 


4.5 (2/7.) 


9 r tntnc 


kpc 






9 


9 


kpc 


L, 


3.4 X 10" 


3.9 X 10" 


3.0 X 10" 


Lo 


Mi 


2.4 X 10^ 


0.5 X 10' 


1.4 X 10' 


Me 


i 


84.5° 


88° 


88° 




Lcloud 


1.0 X 10* 




3.7 X 10** 


Le 


' V 


0.2 


0.4 


1.1 




h./z. 


9.5 


7.5 


7.5 


8.0 (3.6 /im) 


hdlZd 


17 


15 


30 


10 (8.0 yum) 



scale-length, ~ r,runc, an increase in total dust mass (three 
times the smooth model mass) and some 60% of the GMCs to 
be internally heated (equivalent to ~12% Lb„i of the stellar disk). 

The smooth and partially heated clumpy SED models il- 
lustrate how, in the case of a small disk as NGC 4244, high- 
resolution sub-mm data reveals a dust disk with more cold dust 
than a FIR-only model, which is distributed more in clumps and 
throughout the stellar disk (similar scale-length and height as 
the stellar disk). A complete fit of the SED (including future 
Herschel/PACS observations) will be presented in the context of 
our complete low-mass edge-on sample. 

4. Discussion 

Stacked observa tions in the Herschel H-ATLAS survey 
dEales et aUl2010a ) point to specific characteristics {Td ~ 15 K, 
Md ~ 10' Mq) for blue, low-mass spirals (~ lO^M© stellar mass 



T ' I ' I ' I 




X [ (xm ] 

Fig. 3. The Spectral Energy Distribution of NGC 4244. The 
Herschel/SPIRE fluxes from this work (blac k points) ar e show 
together with literature data (red symbols, ' Dale et aO 1200 9t 
MacLachlan et al., 2011). We show the smooth model (blue 
line); the clumpy model with and without embedded stellar 
emission (solid and dashed green line, respectively), and the 
model from MacLachlan et al. (black dotted line). 



iBourne et all 12012'). N GC 4244's stellar and Hi mass are 6 an d 
2.5 X lO'^Mo (Stricklan detaLl l2004t IZschaechner et all 1201 ll) . 
Compared to the H-ATLAS results, our SED dust mass is rea- 
sonable and the gas-to-dust ratio (535 - 170 for the smooth 
and clumped model respectively) is in Une with those found for 
smaller spiral disks (iLerov et al.L l2009h . 
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The central optical depth of t^"' ~ 0.45 - 1.12 is 
similar to the face-on optical depths in occulting galaxies 
( White et al., 2000; Holwerda et al., 2007b) or distant galaxy 
counts faolwerda et all l2005l l2007allcl) . We therefore feel con- 
fident that our SED model accurately reflects the true dust mass 
of NGC 4244. 

The scale-length of the dust in the SED models is consis- 
tently much larger to the stellar one; a factor 1.8 to 4.9, de- 
pending on t he model. For NGC 891, a ratio of 1.5-2 i s typ- 
icaUy found dXilouris et all 119991: iPopescu et al.', '2000", '201lt 



[Driver et al.L 120071: iBianchil 120081; |Bianchi & Xilouris, 2011). 
Similarly, Holwe rda et al. and Holwerda et al. submit- 

ted report a clumpy, extended and flat dust profile in a similar 
mass face-on galaxy occulting a distant bright bulge (2MASX 
J00482 185-2507365, near NGC 253). These flat radial profiles 
of dust suggest radial mixing of the ISM to a similar leve l found 
in the stellar disks (Ide Jong et al.L 120071: iRoskar et al.l l2008i 
Radburn-Smith et al. submitted). 

We find a dust scale- height similar to the stellar one, as do 
iMacLachlan et al.1 ('2011V They attribute this to either the flat- 
tening of dust disks (Dalcanton et al., 2004) or a smaller stellar 
scale-height due to inefficient dynamical heating of th e stellar 
disk in lower mass galaxies (.Bizvaev & Kaisinl |2004|) . To dis- 
tinguish between these two scenarios, we compare the oblate- 
ness (h/z) of the stellar and dust disk (Table [U to those in 
Figure 2 in Dalcanton et al. (2004). We also include the oblate- 
ness of the 3.6 and 8.0 /im Spitzer observations (Holwerda et al 
in prep.). We note that the val ue for z* of 300 pc is the scale- 
height offfiefliin disk found bv lComer6n et al.1 (1201 lb and RGB 
stars dSeth et al.L l2005h . NGC 4244's stellar disk oblateness is 
at the highest, i.e., flattest value for low-mass disks, similar to 
the oblateness of massive disks, arguing against the less effi- 
cient disk heating compared to massive disks. The dusty disk 
oblateness is extremely high, even compared to the PAH emis- 
sion {hjz - 10 at 8 jum), but this is predominantly because of 
the large scale-length of the dust disk. The similar stellar and 
dust scale-height and the best match b y clumpy SED model ar e 
both consistent with the scenario from iDalcanton et"an (|2004|) ; 
in a low-mass disk such as this one, the ISM does not collapse 
into a thin plane but remains fractured and spread throughout the 
height of the stellar disk. 

5. Summary & Conclusions 

Based on new sub-mm observations with the Herschel Space 
Observatory of NGC 4244, and some illustrative models with 
the TRADING code, we find that: 

i. The cold dust is confined to the inner, flat Hi disk (Fig.[T]). 

ii. The dusty ISM is confined to the stellar radius (Fig. |2]i. 

iii. A smooth model under-predicts 70 & 500 /im fluxes (Fig. [3]). 

iv. A clumpy dust disk model with a fraction of the clumps 
heated by embedded sources fits both the general SED (Fig. 
O and the morphology of the disk (Fig.|2]l in sub-mm well. 

vi. The inferred scales of stars and dust are consistent with a 
scenario that includes some radial migration of dust and a 
vertical ISM scale supported by turbulence (Table[T}. 

Future improvements of the NGC 4244 SED model includes 
new FIR (Herschel/PACS) information to constrain the cold dust 
structure and the number of GMCs and the stellar type of their 
heating sources. We plan to put these results in the context of the 
full NHEMESES and HEROES Herschel observations to infer 
if there is indeed a phase-change with mass in the dust structure 
and additionally to benchmark current SED models. 



The authors thank ESA for the Fellowship Program and op- 
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for his SED data and Laura Zschaechner for her Hi map and the 
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